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oo ; Abstract 

^ ' The Skyrme energy-density functional approach has been extended to study the massive heavy- 

■ ion fusion reactions. Based on the potential barrier obtained and the parameterized barrier distri- 

, bution the fusion (capture) excitation functions of a lot of heavy-ion fusion reactions are studied 

' systematically. The average deviations of fusion cross sections at energies near and above the bar- 

\o : 

' riers from experimental data are less than 0.05 for 92% of 76 fusion reactions with Z1Z2 < 1200. 

, For the massive fusion reactions, for example, the ^^'^U-induced reactions and ^^Ca+^'^^Pb the 

O I capture excitation functions have been reproduced remarkable well. The influence of structure 

C ' effects in the reaction partners on the capture cross sections are studied with our parameterized 



- 

X 



barrier distribution. Through comparing the reactions induced by double-magic nucleus ^^Ca and 



^ [ by ^^S and ^^Cl, the 'threshold-like' behavior in the capture excitation function for ^^Ca induced 

reactions is explored and an optimal balance between the capture cross section and the excitation 
energy of the compound nucleus is studied. Finally, the fusion reactions with '^^S, '^''Cl, ^^Ca and 
50Ti bombarding on ^^^Cm, 247,249^]^^ 250,252,254Qf 252,254^3^ ^gj^ ^ ^j^g reactions lead 

to the same compound nucleus with Z = 120 and = 182 are studied further. The calculation 
results for these reactions are useful for searching for the optimal fusion configuration and suitable 
incident energy in the synthesis of superheavy nuclei. 
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I. INTRODUCTION 



It is of great importance to predict fusion cross sections and to analyze reaction mech- 
anism for massive heavy-ion fusion reactions, especially for fusion reactions leading to su- 
perheavy nuclei. In those reactions, the calculation of the capture cross section is of crucial 
importance. It is known that Wong's formula based on one-dimensional barrier penetra- 
tion can describe the fusion excitation function well for light reaction systems, while it fails 
to give satisfying results for heavy reaction systems at energies near and below the barrier. 
For solving this problem, a fusion coupled channel model was proposed, in which the 
macroscopic Woods-Saxon potential together with a microscopic channel coupling concept 
is adopted. By this model fusion excitation functions of some reactions at energies near and 
below the barrier are successfully described. However, it has been found that the parameters 
in the Woods-Saxon potential greatly influence on the results Is] and for heavy systems the 
potential parameters need to be readjusted in order to reproduce experimental data P|. How 
to determine the parameters is still an open problem for predicting fusion cross sections of 
unmeasured reaction systems. Therefore, it is highly requisite to propose a new method for 
systematically describing fusion reactions from light to heavy reaction systems. 

In our previous paper P], we applied the Skyrme energy- density functional for the first 
time to study heavy-ion fusion reactions. The barrier for fusion reaction was calcu- 
lated by the Skyrme energy-density functional together with the semi-classical extended 
Thomas- Fermi method^]. Based on the interaction potential barrier obtained, we pro- 
posed a parametrization of the empirical barrier distribution to take into account the multi- 
dimensional character of the real barrier and then applied it to calculate the fusion excitation 
functions of light and intermediate-heavy fusion reaction systems in terms of the barrier pen- 
etration concept. A large number of measured fusion excitation functions at energies around 
the barriers were reproduced well. Now we try to extend this approach to study very heavy 
fusion reaction systems which may lead to the formation of superheavy nuclei. In these 
cases, the reaction mechanism is very complicated and the capture process is the firstly 
concerned process, which follows by the quasi-fission and fusion, and then the fused system 
further undergoes fusion-fission and evaporation. 

The study of the fusion mechanism (or capture process in very heavy fusion systems), 
especially of the possible enhancement of the fusion (capture) cross section in neutron- 
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rich reactions and also of the suppression of fusion (capture) cross section induced by the 
strong shell effects of projectile or target, is very interesting and essentially importance for 
the synthesis of superheavy nuclei. For fusion reactions induced by double-magic nucleus 
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usion cross sections 







and 4«Ca+y"'y«Zr 



Ca, there exists a puzzle: on one hand, it has been found that the 
at sub-barrier energies are suppressed in fusion reactions ^^Ca+'^^Ca 
0,1^ compared with the ^°Ca+^^Ca and ^°Ca+^°'^^Zr, respectively. On the other hand, the 
experiments of production of superheavy elements Z = 114 and 116 in "hot fusion" reactions 
with ^^Ca bombarding Pu and Cm targets[lo| indicate that the reactions with ^^Ca nuclei, 
indeed, are quite favorable for the synthesis of superheavy nuclei. Therefore, it is worthwhile 
to explore the puzzle concerning the fusion reactions induced by '^^Ca. For this purpose, the 
influence of shell structure, that is, the influence of the Q-value in fusion (capture) process 
on the fusion (capture) cross section is considered in our approach. The choice of an optimal 
reaction combination and a suitable incident energy is always of crucial importance for the 
synthesis of new superheavy nuclei. In order to choose a suitable incident energy, an optimal 
balance between capture cross section and excitation energy of compound nuclei should be 
taken into account. Thus, in this work, a series of fusion reactions induced by ^^Ca, ^^S, 
^^Cl and ^'^Ti are investigated within our approach and the optimal incident energies for the 
reactions are given. 

II. MICROSCOPIC INTERACTION POTENTIAL BARRIER AND 
PARAMETERIZED BARRIER DISTRIBUTION 

In this section, we briefly introduce our approach for calculating the interaction potential 
barrier and fusion (capture) excitation function, a more detailed description can be found 
in ref.jSj. The nucleus-nucleus interaction potentials of fusion systems are calculated within 
the microscopic Skyrme energy-density functional together with the semi-classical extended 
Thomas- Fermi (ETF2) approach (up to second order of h). The interaction potential Vb{R) 
between reaction partners can be written as 

where R is the center-to-center distance between reaction partners, Etot{R) is the total 
energy of the interaction system, Ei and E2 are the energies of the non-interacting projectile 
and target, respectively. The interaction potential Vf,{R) is also called the entrance-channel 



potential in ref. or fusion potential in r ef.ll^. The Etot{R), Ei, E2 are determined by 
the Skyrme energy-density functional P. Illl [isi llJ. Iisj] . 



Etot{R) = j ^[pip(r) + P2p(r - R), pi„(r) + p2n(r - R)] rfr, 
Ei = j H[pip(r),pi„(r)] rfr, 



Eo 



^[P2p(r),p2n(r)] dr. 



(2) 

(3) 
(4) 



Here, pip, p2p, pi„ and p2n are the frozen proton and neutron densities of the projectile 
and target, and the expression of the energy- density functional Ti can be found in refs.0, 
Once the proton and neutron density distributions of the projectile and target are 
determined, the interaction potential Vb(-R) can be calculated from eqs.(l) — (4). 

By density- variational approach and minimizing the total energy of a single nucleus given 
by the Skyrme energy-density functional 7i, the neutron and proton densities of this nucleus 
can be obtained. In this work we take the neutron {i = n) and proton [i = p) density 
distributions of nuclei as spherical symmetric Fermi functions, 

1 



POi 



1 + exp 



■Oi 



{n,p} . 



(5) 



Only two of the three quantities poj , Roi and Oj in this relation, are independent because 
of the conservation of the particle numbers Ni = J pi{r)dr, Ni = {N, Z}. For example, pop 
can be expressed function of Rop and a^, 

-1 



1 + 7r' 



R 



Op 



(6) 



with high accuracy^] when i?op ^ CLp- By using an optimization algorithm, one can obtain 
the minimal energy and the corresponding i?op; CLp, RomO-n for neutron and proton den- 
sities. Then, with the neutron and proton densities of projectile and target obtained we 
can calculate the entrance-channel potential with the same energy-density functional. For 
systematically investigating massive heavy-ion fusion reactions with a simple self-consistent 



17t |. an optimal balance between the ac- 



manner provided by the density functional theory 
curacy and computation cost is adopted in this approach, which is especially valuable for 
theses cases. 
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FIG. 1: (Color online) The entrance-channel potential for reaction ^^Ca+^'^Zr. 

The Skyrme force SkM* jisl is adopted in this work. For a certain reaction system, the 
entrance-channel potential is calculated in a range from R = 7fm to 15/m with a step 
size AR = 0.25/m. Fig.l shows the entrance-channel potential of ^^Ca+^°Zr. The solid 
and crossed curves denote the results of this approach and of the proximity potential [isj]. 
respectively. The results of Skyrme energy-density functional approach are generally close 
to those of proximity potential in the region where the densities of the two nuclei do not 
overlap. The barrier height Bq, radius Ro and the curvature huo near Rq as well as the 
position of fusion pocket Rs can be obtained from the calculations (see Fig.l). Here, the 
curvature hwQ of the barrier is obtained through fitting the entrance-channel potential in 
the region from Rq — 1.25/m to Rq + 1.25/m by an inverted parabola (if Rq — 1.25fm < R^, 
then from R^ to Rq + 1.25/m). 

To overcome the deficiency of one- dimensional barrier penetration model for describing 
sub-barrier fusion of heavy systems, we take into account the multi-dimensional character of 
realistic barrier jlfll due to the coupling to internal degrees of freedom of the binary system. 
We assume that the one-dimensional barrier is replaced by a distribution of fusion barrier 
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D{B). The distribution function D{B) satisfies: 

D{B)dB = 1. (7) 

Motivated by the shape of the barrier distribution extracted from experiments, we consider 
the weighting function to be a superposition of two Gaussian functions Di{B) and D2{B), 
which read 

[B-B,r 



and 



with 



Di{B) 



7 



exp 



-7- 



D2{B) 



exp 



{B - B2) 

i2W2? 



(8) 



(9) 



^i;i = -(5o-5,) 



(10) 



W2 = 2(^0 - Be 



(11) 



Bi = Be + wi, 



(12) 



B2 = Bc + W2. 



(13) 



Here i?o is the height of the barrier (see Fig.l). The B^ = JBq is the effective barrier height 
with a reducing factor / to mimic the lowering barrier effect which is due to the couphng 
to other degrees of freedom, such as dynamical deformation and nucleon transfer, etc. We 
set the reducing factor / = 0.926 in this work, which is the same as in P]. The quantity 7 
in Di{B) is a factor to taken into account the structure effects, which influences the width 
of the distribution Di{B). For the fusion reactions with non-closed-shell nuclei but near 
the /3-stability line we set 7 = 1; for fusion reactions with neutron-shell-closed nuclei or 
neutron-rich nuclei an empirical formula for the 7 values, used in the weighting function 
Di{B) for systems with the same Zi and Z2, was proposed in ref. 0] as 



7 = 1- CoAQ + 0.5(5^' + 



(14) 
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where AQ = Q — Qo denotes the difference between the Q-values of the system under 
considering for complete fusion and that of the reference system. The reference system, in 
general, is chosen to be the reaction system with nuclei along the /3-stability line [^. The 
value of Co is O.SMel/-^ for AQ < and O.lMel/-^ for AQ > 0. The quantities ST^^^"''^^ 
are 1 for neutron closed-shell projectile (target) nucleus and for non-closed cases. 
The fusion excitation function is then given by 

D{B)aJ::\E,,^„B)dB, (15) 



with 



Wong/ IT, D\ ^oRq 1 f 1 , 

(Tfus [Ecm., B) = — In 1 + exp 



■{Ec.m. — B) 



. V...... ~, (16) 

nhjQ 

where £"0.111. denotes the center-of-mass energy, and B, Rq and hwQ are the barrier height, 
radius and curvature, respectively. Using the parameterized barrier distribution functions 
Di{B) and D2{B), we also can obtain the cross sections o"i(-Ec.m.) and cravr(-£'c.m.) by (15) 
with D{B) taken to be Di{B) and D^„{B) = [Di{B) + D2{B)]/2, respectively. Finally, the 
fusion cross section is given by 

cm. 

) = min[(Ti(i?c .111. ) 7 ^avr (^cm.)]. (17) 

The cross section calculated with (17) is referred to as fusion cross section for a light and 
intermediate-heavy system and as capture cross section for a very heavy system at -Ec.m.- 

III. CALCULATED RESULTS FOR FUSION (CAPTURE) EXCITATION 

FUNCTIONS 

In order to extend our approach to study the fusion reactions leading to superheavy nuclei, 
we first check the suitability and reliability of our description of heavy-ion fusion reactions. 
We calculate the fusion excitation functions of 76 fusion reactions with Z1Z2 < 1200 at 
energies near and above the barrier (with 7 = 1) and their average deviations Xiog fi'oin 
experimental data defined as 



A. log 

n=l 



^ m 

- J2 [log(ai,(K)) - log(ae,.,(K))]' . (18) 



Here m denotes the number of energy-points of experimental data, and ath{En) and <Jexp{En) 
are the calculated and experimental fusion cross sections at the center-of-mass energy E^, 
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FIG. 2: (Color online) The average deviations for a total of 76 fusion reactions with Z1Z2 < 



1200. The solid circles and crosses denote t^ 
potential with fixed potential parameters 



le results of our approach and those with a Woods-Saxon 
respectively. In the calculations of fusion cross sections 
at energies near and above the barrier with the Woods-Saxon potential, the code CCFULLy] is 
used without taking into account the excitation and deformation of the reaction partners. 

{En > Bq), respectively. Fig.2 shows the results for X\og which the solid circles and 
crosses denote the calculated results from this approach and those from ref. Q], respectively. 
Applying the approach of ref.P| there are 43% systems in 76 fusion reactions in which the 
average deviations Xiog of calculated fusion cross sections from the experimental data are less 
than 0.05, but for reactions with Z1Z2 > 640 the results are not satisfying very well. With our 
approach, the average deviations of 92% systems in Xiog l^ss than 0.05, which indicates 
that this approach is successful for describing fusion cross sections of heavy-ion reactions 
at energies near and above the barrier from light to intermediate-heavy fusion systems. In 
Fig. 3 we show three examples of fusion excitation functions for the reactions ^^0+^^'*Sm[ 



160+922r 



211 and ^^Ni+^^Zrl2a|, in which the solid and dashed curves present the results of 

Q 

our approach and of ref. [2|, respectively, the squares denote the experimental data. From 
this figure we can see that our approach gives quite reasonable description for all selected 
fusion reactions with the Z1Z2 up to 1120 at energies near and above the barrier. 
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FIG. 3: (Color online) Fusion excitation functions for ^^O+^^^Sm, and ^'^m+^'^Zr. The 

squares and solid curves denote the experimental data and the results of this work, respectively. 
The dashed curves denote the results of the approach with the Woods-Saxon potential with fixed 
potential parameters [3]. 

For more massive fusion reactions leading to superheavy nuclei, the quasi-fission process 
occurs and therefore, the capture cross sections are larger than the corresponding fusion 
cross sections. In ref.j2^ the fission and quasi-fission process in ^^^U-induced reactions were 
studied. Fig. 4 shows the results in which the solid and open circles denote the measured cross 
sections for the fission-like process and for complete fusion followed by fission, respectively. 
The solid curves give the calculated results of our approach with 7 = 1. From this figure 
one can see that the calculated capture excitation functions of the reactions ^^^U+^^Mg, 
^^^11+^'^ Al, ^^^U+^^S, and ^^^U+^^Cl are quite close to the measured fission-like cross sec- 
tions. It implies that our approach can describe the massive fusion reactions between nuclei 
with neutron open shells but near the /5-stability line. 

For the very massive fusion reaction between double-magic nuclei '^^Ca and ^°^Pb, the 
infiuence of the shell effects is very significant. So careful consideration of the 7 value is 



required at sub-barrier energies. Fig. 5 shows the calcu 



^^Ca+^°^Pb and the experimental data of refs 



.y and Se- 



ated capture excitation function of 



The dashed curve presents the 



results with 7 = 1, that is, no neutron-shell-closure effect is considered. The solid curve is 
calculated with 7 = 9.5 according to eq.(14), in which the closed shell effect is considered. 
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FIG. 4: (Color online) Capture cross sections of ^^^U+^^Mg,^''Al,^^S,^^Cl. The solid and open 
circles denote the measured cross sections for a fission-like process and for complete fusion followed 
by fission, respecti vely . The solid curves are the results from our approach with 7 = 1. The stars 



are taken from ref. 



24|. 



We find that for energies below the barrier the experimental data can only be described 
with 7 = 9.5 and the calculations with 7 = 1 are over-predicted. From this analysis, one 
learns that the measured capture cross sections of ^^Ca+^'^'^Pb at sub-barrier energies are 
obviously suppressed, which may arise from the suppression of the nucleon transfer between 
reaction partners due to the strong closed shell effects, which will be further studied in the 
following section. 

IV. OPTIMAL BALANCE BETWEEN CAPTURE CROSS SECTION AND 
EXCITATION ENERGY OF COMPOUND NUCLEI 

It is very important to find a favorable combination of projectile and target and a suitable 
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FIG. 5: (Color online) Capture cross sections of ^^Ca+^'^^Pb. The solid and open circles denote 



and ref. 



respectively. The dashed and 



the measured capture-fission cross sections from ref. 
solid curves present the calculated results with 7 = 1.0 and 9.5 obtained by eq.(14), respectively. 
The dash-dotted line indicates the energy corresponding to the height of the barrier. 

incident energy for synthesis of superheavy nuclei. In this section we study very massive 
fusion reactions and search for an optimal balance between the capture cross section in the 
entrance channel and the excitation energy of the compound nuclei. For searching a fusion 
system with large capture cross sections, we carried out a series of calculations for fusion 
reactions induced by ^^'^^S, 35,37qj ^^^^ ^^Ca projectiles. For example, Fig. 6 shows the capture 
excitation functions of the reactions ^^S-|-^^^Cf and ^^Cl-|-^^^Es. The solid curves present the 
results with 7 = 1 (without considering structure effects in the entrance channel), and the 
dashed curves are for the results with the 7 obtained from (14) i.e. 7 = 0.5 for ^^S-|-^^^Cf 
and 7 = 0.6 for ^^Cl-|-^^^Es. The enhancement of capture cross sections in the sub-barrier 
energy region with the 7 < 1 is caused by the effect of excess of neutrons in reaction systems. 
So from the point of view of increasing the capture cross sections, it is more favorable to 
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FIG. 6: (Color online) Capture excitation functions for fusion systems ^^S+^^^Cf and ^^Cl+^^^Es. 
The dash-dotted lines indicate the corresponding barriers. The solid and dashed curves denote the 
results with 7 = 1 and with 7 value obtained with eq.(14), respectively. 

select the reaction systems with 7 < 1. However, the amount of the excitation energy of 
the formed compound nucleus is essentially important for the survival probability. The 
smaller the excitation energy is, the larger the surviving probability is. Thus, seeking an 
optimal balance between the capture cross section and the excitation energy of the compound 
nucleus becomes very important for synthesis of superheavy nuclei . For choosing the fused 
nuclei with an excitation energy as low as possible, the fusion reactions with double-magic 
nuclei ^^Ca are considered to be good candidates because of the low Q- values for those 
fusion reactions. As an example, let us investigate reaction ^^Ca+^^^Cm. For this reaction 
the 7 value is equal to 10.8 calculated with eq.(14). Fig. 7 shows the capture excitation 
function for this reaction, in which the solid and dashed curves denote the results for the 
cases of 7 = 1 and 7=10.8, respectively. From this figure one finds that for fusion reactions 
induced by double-magic nuclei ''^Ca the capture cross sections at sub-barrier energies are 
suppressed compared with reactions with shell-open nuclei but near the /^-stability line. 
However, if we suitably choose an incident energy, for example, as indicated by the arrow 
in Fig. 7, the capture cross section of the reaction ^^Ca+^^^Cm is not suppressed so much 
(still reaches several tens of milli-barns) and the excitation energy of the compound nuclei 
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FIG. 7: (Color online) Capture excitation functions for ^^Ca+^^^Cm. The solid and dashed curves 
represent the results with the 7 =1 and with 7 =10.8 obtained from eq.(14). The arrow indicates 
the incident energy at which the corresponding excitation energy of the formed compound nucleus 
is E^j^ = 3lMeV. 

is only E^j^ = SlMeV. Such an incident energy was already used in the experiment of 
ref.|lO|. Now let us make a comparison between the reaction ^'^Ca+^'^^Cm and the reactions 
32g_|_254Qf g^j^^ ^^Cl+^^^Es. For the system '^^Ca+^^'^Cm, the capture cross section is of about 
80mb and the excitation energy is about 31MeV if the incident energy is taken to be about 
198MeV. While, for the systems ^^S+^^^'Cf and ^^Cl+^^^Es, if the same excitation energy is 
required the incident energies must be as low as about 150MeV and IGOMeV, respectively, 
since the Q- values of these two fusion reactions are much higher compared with ^^Ca induced 
reactions. At these incident energies the capture cross sections for these two reactions are as 
small as those less than O.lmb according to this model calculations. From above analysis we 
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FIG. 8: (Color online) Fusion barrier distribution for ^^O+^'^^Pb. The distribution is evaluated 
with AEc.m. = 2.5MeV . The solid squares and solid curve present the experimental data and the 
results from our calculations, respectively. 

can conclude that the fusion reaction ^^Ca+^^^Cm seems to be more favorable compared to 
32g_|_254Q£ ^^Cl+^^'^Es if a suitable incident energy is chosen, as far as both the capture 
cross section and the excitation energy of the compound nuclei are concerned. 

Now let us discuss how to choose suitable incident energy. We notice that the capture 
excitation function for reactions induced by double-magic ^^Ca goes very sharply down at 
sub-barrier energies due to strong closed shell effects, as shown by the dashed curve of Fig. 7. 
It seems to us that there exists a 'threshold-like' behavior, which is important for choosing 
the incident energies. This 'threshold-like' behavior of excitation function of capture cross 
sections is closely related to the shape of the barrier distribution. In our previous paper 



a 



a number of barrier distributions were calculated according to expressions (8_ 



•(13). As 
3l| in Fig.8. 



example, here we show the calculated fusion barrier distribution for ^^0+^°'^Pb 
The agreement of the calculated barrier distribution with experimental data tells us that 
our approach about the parameterized barrier distribution is quite reasonable. The effective 
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weighting function De«{B) is defined as 

f Di(B) : B <B^ 
Des{B)=l ' (19) 

[Davr(5) : B>B, 

(with / Dcs{B) dB ^ 1 and J D^^^{B) dB = 1, see ref.[5]). The B^ denotes the position 
of the left crossing point between Di{B) and Davv{B). The function D^^iB) can describe 
the fusion excitation function reasonably well. Fig. 9 shows the capture excitation function 
(Fig. 9(a)) and the effective weighting function D^s (Fig. 9(b)) for the reaction ^^Ca+^^^Pu. 
The dotted vertical line denotes the barrier height Bq, and the short dashed vertical line 
indicates the energy at the peak of Des which we call the most probable barrier height -Bm.p.- 
From the dashed curve of Fig. 9 (a) one can see that the capture cross section goes down very 
sharply when the incident energy is lower than -Bm.p.- This is because the decreasing slope 
of the left side of the weighting function Dcs is very steep due to strong closed shell effects 
(7 = 11.0). In fact, one can find that the left side of the barrier distribution D^s^B) is 
given by Di{B) (see expression (19)), which becomes a 5-function when 7 00. For the 
system with 7 much larger than 1 the effective barrier Dgfr will have the similar behavior 
as is shown in Fig.9(b). Thus, the most probable barrier energy B^ p can be considered as 
the incident energy 'threshold', and for massive fusion reactions with 7 much larger than 
1 leading to superheavy nuclei such as ^*Ca induced reactions, the suitable incident energy 
should be chosen in the region Ec.m. > -Bm.p.. The barrier distribution for this case shown in 
Fig. 9(b) looks like a 5-function with a long tail in the high energy side. It seems to be that 
the Wong's formula with barrier height being the -Bm.p. should work without introducing 
the 7. But the results calculated with Wong's formula and expression (17) are different 
especially at sub-barrier energies, as shown in Fig. 9 (a) (comparing the dot-dashed curve 
and the dashed curve). It seems to us that with this 7 value like 7 = 11.0 the behavior of 
-Doff is still different from a 5-function and the parameter 7 still plays a role. 

We find that the incident energies adopted in the experiments successfully producing 



superheavy nuclei in recent years jSc 



y,Q,Ql 



for some reactions induced by ^®Ca are 



very close the most probable barrier energies -Bm.p.. Table 1 gives the comparison of the 



calculated most probab 



used in recent years 



carrier energies -Bm.p. with experimental incident energies E[ 



)able barrier t 

3, 24, HQ 



exp 



for some reactions induced by ^^Ca leading to producing 



superheavy nuclei. The barrier hight Bq, position Rq of the barrier, curvature at the top of 



15 




FIG. 9: (Color online) (a) Capture excitation function and (b) effective weighting function for the 
reaction ^^Ca+^^^Pu. In (a) the solid and dashed curves show the results with 7 = 1 and with 7 
obtained by eq.(14), respectively. The dot-dashed curve denotes the results from Wong's formula 
with B = -Bm.p.- The results in (b) are obtained by setting 7 = 11.0. 
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FIG. 10: (Color online) Capture excitation functions for the systems (a) ^^S+^^^Es and (b) 

48Ca+254Es. 

the barrier expressed by huo, factor 7 are also listed. In addition, we list the mean value 
-Bmean of the barrier height defined as 

_ jBD^^jB) dB 
~~ jD,siB)dB ■ ^ ^ 

The -Bmean IS, in general, larger than the -Bm.p. since the slope of the left side of the weighting 
function Des is very steep. From the table one can find that for all listed reactions the 
energies E!^[^ are higher than the calculated most probable barrier energies -Bm.p., which 
supports our ideas about how to choose the favorable incident energy. Further, we find the 
experimental evaporation-residue excitation functions of the fusion reactions listed in Table 
1 are peaked at the energies ranging from Smean to Bq in most cases, which implies that the 
energy i?mcan niay be more suitable to be chosen as the incident beam energy in the fusion 
reactions with 7 much larger than 1 for producing superheavy nuclei. 

In addition to the reactions induced by ^'^Ca leading to superheavy nuclei, reactions with 
36S, 37 CI, 48Ca and ^OTi bombarding on ^^^Cm, 247,24931^^ 250,252,254^^ ^^^^ 252,254^5 ^^^^ ^^^^ 

studied and all relevant parameters for the entrance-channel capture barriers for those fusion 
reactions are listed in Table 2. The table gives the Q-value for the reactions, the barrier 
height Bq, position Rq of the barrier, curvature at the top of the barrier expressed by hwQ, 
factor 7 of structure effects, the mean value i?mean of the barrier, the most probable barrier 
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TABLE I: The entrance-channel capture barriers of fusion reactions with Ca nuclei. 



Reaction B^{MeV) Ro{fm) HwoiMeV) 7 B^,^^{MeV) B^,^,{MeV) E'^l{MeV) 



48ca+207pb [33j 


183.37 


12.0 


4.44 


8.9 


176.28 


173.18 


173.3 


48ca+208pb [33j 


183.17 


12.0 


4.43 


9.5 


176.10 


173.03 


173.5 


48ca+238u [34] 


200.82 


12.25 


4.19 


10.7 


193.09 


189.71 


191.1 


48Ca+242pu [3^ 


204.78 


12.25 


3.90 


11.6 


196.91 


193.44 


196 


48Ca+244pu ^ 


204.31 


12.25 


3.99 


11.0 


196.44 


192.99 


193.3 


48Ca+243Am [36j 


206.87 


12.25 


3.87 


9.8 


198.89 


195.39 


207.1 


48ca+245cm [35j 


208.80 


12.25 


3.88 


11.7 


200.77 


197.22 


203 


48ca+248cm [34] 


208.25 


12.25 


3.89 


10.8 


200.23 


196.71 


198.6 



energy -Bm.p., the excitation energy of compound nucleus E^,^ when -Ec.m. = -Bmean, and the 
depth of capture pocket Bq—Bs (or called quasi- fission barrier height j37|, here Bg denotes the 
value at the bottom of the pocket, see Fig.l). Comparing the data from different reactions 
one can find that the reactions with ^''Cl induce relatively higher excitation energies E'^^ and 
those with ^^Ca and ^°Ti produce relatively lower excitation energies when -Ecm. = -Bmean- 
So ^^Ca and ^°Ti induced reactions can be considered as good candidates of cold fusion 
reaction for producing superheavy nuclei from the point of a low excitation energy of the 
compound nuclei. Here we have not studied the orientation effect of deformed target, which 
has significant effect on fusion barrier height and the compactness of the fusion reactions. 
Recently, compactness of the ^^Ca induced hot fusion reactions was studied in which it was 
shown that ^^Ca induced reactions on various actinides were the best cold fusion reactions 
with optimum orientations of the hot fusion process 33]. By comparing the depths of the 
capture pockets for different reactions we find that the depth decreases with increase of 
the proton number of the projectile nuclei. We know that the shallower the pocket is, the 
stronger the quasi-fission is. So the projectile ^^S inducing capture reactions is more favorable 
for the small quasi-fission probabilities of those reactions. By using this table we can easily 
calculate the capture cross sections by eqs.(15) — (17) for all reactions listed. Fig.lO shows the 
calculated capture excitation functions for the systems ^^S+^^^Es and ^^Ca+^^^Es with our 
approach by using the data from Table 2. In addition, the entrance-channel capture barriers 
of the reactions ^^Ni+^^su, ^^Fe+^^^Pu, ^^Ci+'^^^Cm and 50Ti+252Cf which lead to the same 
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compound nucleus with Z = 120 and = 182 are calculated and listed in Table 3. Table 2 
and Table 3 provide us with very useful information for choosing an optimal combination of 
projectile and target and suitable incident beam energies for producing superheavy nuclei 
for unmeasured massive fusion reactions. 

V. CONCLUSION AND DISCUSSION 

In this work, the Skyrme energy-density functional approach has been applied to study 
massive heavy-ion fusion reactions, especially those leading to superheavy nuclei. Based on 
the barriers calculated with the Skyrme energy-density functional, we propose the parame- 
terized barrier distributions to effectively taken into account the multi-dimensional character 
of the realistic barrier. A large number of heavy-ion fusion reactions have been studied sys- 
tematically. The average deviations of fusion cross sections at energies near and above 
the barriers from experimental data are less than 0.05 for 92% of 76 fusion reactions with 
Z1Z2 < 1200. Massive fusion reactions, for example, the ^^^U-induced reactions and the 
'^^Ca+^^'^Pb have been studied and their capture excitation functions have been reproduced 
well. The influence of the structure effects in the reaction partners on the capture cross 
sections are studied by using parameter 7 in our model. To search the most favorable con- 
dition for the synthesis of superheavy nuclei, the optimal balance between the capture cross 
section and the excitation energy of the formed compound nuclei is studied by comparing 
the fusion reactions induced by the double-magic nucleus '^^Ca and by ^^S and ^^Cl. Based 
on this study, the 'threshold-like' behavior of excitation function of capture cross sections 
with respect to incident beam energy has been explored and possible values of this 'thresh- 
old' for reactions mainly induced by ^^Ca are given. Finally, we have further studied the 
capture reactions leading to superheavy nuclei such as ^^S, ^^Cl, ''^Ca and ^°Ti bombard- 
ing on 248Cm, 247,24931^^ 250,252,254Qf 252,254^3^ ^^^^ ^J^g rcactloUS 64Ni+238U, 

^^Fe+2'^^Pu, ^^Cr+^^^Cm and ^^Ti+^^^Cf which lead to the same compound nucleus with 
Z = 120 and N = 182. The relevant parameters for calculating the capture cross section of 
these reactions have been provided which are helpful for the study of unmeasured massive 
fusion reactions. Especially, we predicted optimal fusion configuration and suitable incident 
beam energies for the synthesis of superheavy nuclei. 

We notice that the deformation and orientation of colliding nuclei have a very significant 



role on fusion reactions. In 



the effect of deformation and orientation on the barrier 
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hight and the compactness of fusion reactions were investigated systematically. However, 
this kind study is beyond the scope of present work. We have only made preliminary 
calculations of the potential barrier for ^^Ca+^^^Cm with the deformation and orientation 
of ^^^Cm taken into account in the entrance channel. For this reaction the lowest barrier 
is obtained for the orientation = 0°, i.e. when ^^Ca touches the tip of deformed ^^^Cm 
target, while the highest barrier is obtained for © = 90°, when ^®Ca touches the side. The 
lowest barrier obtained for © = ° is a little bit lower than the most probable barrier height 
-Bm.p. of this reaction given in Table 1 and the barrier distribution due to the orientation of 
^^^Cm is close to the effective weighting function DesiB) which is for describing the capture 
process of the reaction if assuming the orientation probability decreases gradually from 0° 
to 90°. So the deformation effects seem to be partly involved in the parameterized barrier 
distribution functions. The study on this aspect is in progress. 
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TABLE II: The entrance-channel capture barriers for fusion reactions with S, CI, Ca and 
50Ti bombarding on ^^SCm, 247,2493^, 250,252,254cf ^nd ^^^'^^^Es. 



Reaction Q{MeV) Bo{MeV) Ro{fm) ^o(MeF) 7 5mean(MeF) Bra.p.{MeV) E*^^{MeV) Bq-Bs 



36(3 1 248/^™ 
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4.o4 
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TO 
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360 1 247t31 
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OT 
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CI 
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48Ca+247Bk 


-171.71 


210.80 


12.25 


3.95 


9.6 


202.67 


199.11 


30.95 


5.27 


48Ca+249Bk 


-170.76 


210.46 


12.25 


3.83 


9.1 


202.33 


198.76 


31.57 


5.30 


48Ca+250cf 


-174.53 


212.56 


12.25 


3.66 


11.4 


204.39 


200.81 


29.86 


5.11 


48Ca+252cf 


-173.77 


212.10 


12.5 


4.41 


11.0 


203.94 


200.37 


30.17 


5.17 


48Ca+254cf 


-173.28 


211.63 


12.5 


4.38 


10.7 


203.48 


199.92 


30.20 


5.25 


48Ca+252Es 


-177.43 


214.29 


12.5 


4.43 


10.6 


206.04 


202.39 


28.61 


4.98 


48Ca+254Es 


-176.97 


213.94 


12.5 


4.39 


10.3 


205.70 


202.13 


28.73 


5.05 
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Reaction Q{MeV) Bo{MeV) Roifm) hwo{MeV) 7 B^^^n{MeV) B^,pXMeV) E^^{MeV) Bq-B, 



50Ti+248Cj^ 


-185.52 


229.00 


12.25 


3.75 


3.9 


219.88 


216.39 


34.36 


4.41 


50Ti+247Bk 


-191.42 


231.85 


12.25 


3.80 


4.1 


222.63 


219.00 


31.21 


4.14 


50Ti+249Bk 


-189.78 


231.45 


12.25 


3.67 


3.3 


222.16 


218.71 


32.38 


4.18 


Tl+ Cl 


-194.40 


233.79 


12.25 


3.64 


4.9 


224.56 


220.84 


30.16 


3.98 


50Ti+252cf 


-193.02 


233.23 


12.25 


3.53 


4.2 


223.97 


220.33 


30.95 


3.96 


50Ti+254cf 


-191.92 


232.67 


12.5 


4.38 


3.6 


223.38 


219.81 


31.46 


3.97 


50Ti+252Es 


-197.90 


235.72 


12.25 


3.52 


4.3 


226.37 


222.67 


28.47 


3.71 




-190.81 


235.21 


12.5 


1.39 


3.8 


225.86 


222.21 


29.05 


3.69 



TABLE III: The same as Table 2, but for reactions ^^Ni-F^ssu, ^^Fe-F^^^Pu, ^^Cr-F^^SCm and 

50Ti+252cf. 



Reaction Q(MeF) Bo{MeV) Roifm) hLOo{MeV) 7 5mcan(MeF) B„,,pXMeV) E^j^{MeV) Bq-Bs 



64Ni+238u 


-237.41 


276.01 


12.5 


4.61 


7.1 


265.26 


260.73 


27.85 


1.79 


58Fe+244p^j 


-219.97 


262.88 


12.25 


3.93 


1.0 


251.56 


248.80 


31.60 


2.56 


54cr+248cm 


-207.16 


248.52 


12.25 


4.20 


3.0 


238.51 


234.86 


31.35 


3.26 


50Ti+252cf 


-193.02 


233.23 


12.25 


3.53 


4.2 


223.97 


220.33 


30.95 


3.96 
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